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Abstract 

The paper presents a multi-body analysis of the 
1/5 scale wind tunnel model of the V-22 tiltro- 
tor, the Wing and Rotor Aeroelastic Testing 
System (WRATS), currently tested at NASA 
Langley Research Center. An original multi- 
body formulation has been developed at the Di- 
partimento di Ingegneria Aerospaziale of the Po- 
litecnico di Milano, Italy. It is based on the di- 
rect writing of the equilibrium equations of in- 
dependent rigid bodies, connected by kinematic 
constraints that result in the addition of alge- 
braic constraint equations, and by dynamic con- 
straints, that directly contribute to the equilib- 
rium equations. The formulation has been ex- 
tended to the simultaneous solution of interdis- 
ciplinary problems by modeling electric and hy- 
draulic networks, for aeroservoelastic problems. 
The code has been tailored to the modeling of 
rotorcrafts while preserving a complete general- 
ity. A family of aerodynamic elements has been 
introduced to model high aspect aerodynamic 
surfaces, based on the strip theory, with quasi- 
steady aerodynamic coefficients, compressibil- 
ity, post-stall interpolation of experimental data, 
dynamic stall modeling, and radial flow drag. 
Different models for the induced velocity of the 
rotor can be used, from uniform velocity to dy- 
namic inflow. A complete dynamic and aeroe- 
lastic analysis of the model of the V-22 tiltrotor 
has been performed, to assess the validity of the 
formulation and to exploit the unique features of 
multi-body analysis with respect to conventional 
comprehensive rotorcraft codes; These are the 
ability to model the exact kinematics of mechan- 
ical systems, and the possibility to simulate un- 


usual maneuvers and unusual flight conditions, 
that are particular to the tiltrotor, e.g. the con- 
version maneuvre. A complete modal validation 
of the analytical model has been performed, to 
assess the ability to reproduce the correct dy- 
namics of the system with a relatively coarse 
beam model of the semispan wing, pylon and ro- 
tor. Particular care has been used to model the 
kinematics of the gimbal joint, that characterizes 
the rotor hub, and of the control system, con- 
sisting in the entire swashplate mechanism. The 
kinematics of the fixed and the rotating plates 
have been modeled, with variable length con- 
trol links used to input the controls, the rotat- 
ing flexible links, the pitch horns and the pitch 
bearings. The investigations took advantage of 
concurring wind tunnel test runs, that were per- 
formed in August 1998, and allowed the acquisi- 
tion of data specific to the multi-body analysis. 


Introduction 

Traditional analysis of rotorcraft is based on 
well established methods, usually implemented 
in dedicated analysis tools. Transfer Matrix, or 
Myklestad, and FEA are used to determine the 
structural properties of flexible blades, such as 
frequencies and mode shapes. Dedicated formu- 
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lations are used to study the basic behavior of ro- 
tors in trimmed hover and forward flight condi- 
tions. Examples are PASTA [13], DYN4-DYN5 
[18]. Usually these formulations rely on basic 
assumptions that simplify the problem, thus im- 
proving the efficiency of the analysis and re- 
ducing the computational cost. On the other 
hand, usually, only a limited set of problems, 
with little flexibility in the choice of the configu- 
ration, can be effectively handled. Comprehen- 
sive codes give an important rotorcraft analysis 
capability. For years they have been the best 
trade-off between both the generality of the for- 
mulation and of the range of application, and 
reasonable simplifications. A relatively detailed 
structural description is usually allowed; the de- 
grees of freedom are successively condensated 
by modal reduction, and the reduced model is 
used for subsequent aeroelastic analyses, ranging 
from trimmed solutions to aeroelastic stability. 
Examples are CAMRAD [10], [11], UMARC [9], 
In some cases, a sophisticated aerodynamic anal- 
ysis is used to determine the wake of the rotor, 
as in CAMRAD, UMARC, CAMRAD II [12], 
A general purpose modeling code represents an 
analysis tool allowing the handling of a wide 
spectrum of problems with as little limitations 
as possible; the multi-body interdisciplinar ap- 
proach is a clear example. A multi-body model 
has unique features with respect to more con- 
ventional approaches: there are no kinematic as- 
sumptions or simplifications. The kinematic be- 
havior of a mechanism can be modeled to the de- 
sired degree of refinement, with exact kinematic 
relationships between bodies; elastic bodies can 
be modeled with a degree of refinement compa- 
rable to that of a nonlinear FEA; the designer 
is left a complete freedom on the description of 
the system, which is built from scratch compo- 
nent by component. This approach is likely to 
need more computer power than that required 
by specialized and simplified approaches, but 
pays back in terms of efficiency since it allows 
the designer to avoid risky physical oversimpli- 
fications along with the greater modeling confi- 
dence allowed by using a single, general purpose, 
and well proven modeling tool. Moreover, with 
the computer power nowadays available, even 
the most complex models are likely to require a 
turnaround time that is compatible with an ex- 
tensive set of parametric analyses and, in a short 
time, even with a complete system optimization. 
The technology of multi-body is not mature yet, 
with particular regard to the aeronautical and 


specifically the aeroservoelastic field. Current 
commercial general purpose multi-body analysis 
codes, e.g. DADS [8], MECANO [3],' ADAMS 
and others [20], still pose some limitations to 
the modeling of rotorcrafts, mainly due to in- 
sufficient aerodynamics, insufficient description 
of flexible bodies, and in some cases to limita- 
tions in the integration algorithms when applied 
to large finite rotations of the order of some rev- 
olutions [14], A special case is represented by 
CAMRAD II, which should be considered some- 
thing more than a comprehensive code, but not 
yet a general purpose, global modeling multi- 
body code. In the research field, a promising 
multi-body code, successfully applied in rotor- 
craft analyses, is DYMORE [1], A brief de- 
scription of the multi-body formulation applied 
in this paper is presented first, followed by a 
description of the tiltrotor model form a multi- 
body standpoint. Finally, the results of the anal- 
ysis are discussed. 

Multi-Body Formulation 

Dynamics. The multi-body problem is formu- 
lated by directly writing the equilibrium equa- 
tions of each body. Constraints are imposed 
by adding constraint equations, resulting in un- 
known “reacting forces” as unknowns, in a La- 
grangian Multipliers style. The reaction un- 
kowns are directly the reaction forces; this over- 
comes the need of post processing the multipliers 
to determine the reactions. The dynamics prob- 
lem has been written as a first order differential 
system of equations; the equilibrium equations 
of a body are: 

f /j — F (.t, ;r, R, u, t) = 0 

\ 7 — (oj x S) x x — M ( x , x, R, u>,t) = 0 

where 0, 7 are the momenta, x is the position 
of the node, R is the rotation matrix that de- 
scribes the rigid rotation from the local to the 
global frame, expressed in terms of the Gibbs- 
Rodriguez rotation parameters <7, to is the angu- 
lar velocity of the node, related to the rotation 
matrix by wx = RR T where operator x rep- 
resents the cross product; to is related to the 
derivatives of the parameters by the expression 
to = Gg. The inertial forces and moments bal- 
ance in a d’Alembert sense the forces F and mo- 
ments M, which may depend on the configura- 
tion, e.g. the elastic forces, as well as on other 
parameters, like the time t. The term (w x S) xx 



in the moment equation is due to the motion of 
the pole the moments are referred to. The defi- 
nitions of the momenta are: 

J P = rnx — S x to 

| 7 = 5xi + Ju> 

The mass of the body is m; the inertial proper- 
ties S, J represent the first and second order in- 
ertia moments of the body, referred to the global 
frame; their transformation from the local to the 
global frame is S = RS, J = RJR T . Kinematic 
constraints are added as constraint equations: 

$ (x, x, R,w , . . . , t) =0 

The unknown constraint reactions Vf, Vm 
contribute to the equilibrium equations as 
Lagrange multipliers. Equations $ may repre- 
sent holonomic (algebraic) and non-holonomic 
(differential) constraints. The resulting system 
of equations, made of the equilibrium equa- 
tions, the definitions of the momenta and the 
constraint equations, is known as Lagrangian of 
the first kind, and it is Differential Algebraic 
(DAE) of index three [2]. The system is solved 
in the unknowns X, g , ft, 7, Vf, Vm without any 
further substitution. Control and Servosystem 
equations are added as generic differential equa- 
tions, with generic scalar unknowns, they are 
implemented in the so called General elements, 
or Genel. 

Kinematics. The formulation of the kinemat- 
ics of finite rotations is fundamental in a multi- 
body implementation. The orientation of a local 
frame is described by an orthogonal matrix R 
that maps vectors from the local to the global 
frame. The parametrization of large rigid ro- 
tations requires at least three unknowns, but 
four parameters are needed to avoid singulari- 
ties that arise when the orientation of the rota- 
tion is undefined. This problem has been pre- 
vented by considering incremental unknown ro- 
tations. The current orientation of the generic 
reference frame R is accounted for by a constant 
rotation matrix R r , multiplied by the unknown 
incremental rotation represented by matrix i?A , 
assumed to be small enough to avoid any sin- 
gularity. This assumpion is reasonable since a 
limited rotation at each step is required for ac- 
curacy considerations. Matrix R r is updated 
at each step. The Gibbs-Rodriguez parameters 
g = 2 tan -1 (<p/2) represent a very efficient fi- 
nite rotations parametrization from a computa- 
tional standpoint, since they are entirely alge- 


braic. The linearized expressions of the rota- 
tional entities concide with those of the rota- 
tional vector ( p . The rotation matrix R is: 

R = I + , , 4 T (ffx +^ff X gx ) 

4 + g J g \ ^ J 

while matrix G is: 

G = Jr?y( 7+ ?' x ) 


Integrator. Time integration is performed by 
an implicit, A/L- stable, second order accurate 
predictor-corrector integrator. The basic formu- 
las are: 

12 

Vk = — (Vk-i ~ Vk- 2) + + 5y k -2 

Vk = {l-a)y k -i+ay k -2+h(d+^jy k 

+ h a + - - 2 y k -i + h a + d'j y k -2 

for the predictor, consisting in the cubic extrap- 
olation of the derivatives based on the states and 
their derivatives at the two preceding steps, and 
the prediction of the state based on an extrapo- 
lation, whose coefficients ensure second order ac- 
curacy, with user-defined control of the algorith- 
mic damping; h is the time step. The formulas 
have been generalized to a variable step predic- 
tor. The coefficients a and <5 can be expressed in 
terms of the desired asymptotic spectral radius 
Poo , under the assumption of real and coincident 
asymptotic roots: 

4 f4 - (! -Poo) 2 x (1 - Poo) 2 

4 — (1 — Poo)” 2(4— (1— Poo) 2 ) 

For poo = 1 the method is very similar to 
the Crank-Nicholson rule (no dissipation), 
though using two steps, while for poo = 0 
the method coincides with the well known 
Backward Difference Formulas (BDF) [2]. The 
correction is performed by a complete/modified 
Newton-Raphson iteration. 

Rotations Updating. The predicted config- 
uration of the system is used as reference for 
the correction in a way called updated- updated. 
Only the rotation related to the correction is un- 
known; as a consequence, it is expected to be 



really small, provided the prediction is accurate. 
The angular velocity becomes: 

to = G aPa + R,AtO r 

where subscripts (-) A underline that the rota- 
tion parameters and their derivatives are re- 
ferred only to the correction rotation. As a con- 
sequence, when the prediction is accurate, the 
terms involved in the linearisations can be ap- 
proximated as: Ra — I, Ga — I, u — Qa + to r . 
The linearisations become: 

A R = Acja x R r AGa = 0 

Aw = Ac/a + Aqa x to r 

This greatly simplifies the writing of the 
Jacobian matrix, with consequent savings in 
computational time; the accuracy is preserved 
by consistently calculating the residual. 

Beams. Beams are the main elastic element 
of the presented formulation. The strains are 
defined as the difference between the current and 
the initial derivatives of the reference line p (£) 
that describes the position of the beam. The 
strains, in the material frame, are: 

£ = r t p'-p ' 0 

where the position p refers to the current frame, 
while po refers to the initial configuration of the 
beam in the material frame, and the prime (•/ 
performs a spatial derivative with respect to an 
abscissa £ along the reference line. The geomet- 
ric curvature p is defined as the spatial derivative 
of the reference frame of the beam section: 

px = R t R' 

The difference between the current and the ini- 
tial, or imposed, curvature po, represents the 
elastic curvature k, again in the material frame: 

kx = R t R' - R^Rq = p x -pox 

When incremental rotations are considered, the 
elastic curvature becomes: 

k = R T Gg' + k r 

The above mentioned simplifications descending 
from the updated-updated rotations approach 
also apply to the beam kinematics. An origi- 
nal finite volume approach is used to formulate 
the beam element. Finite volumes applied to 


the equilibrium of a beam can be physically in- 
terpreted as the direct balance of the forces and 
the moments that act on a finite piece of beam, 
including the internal forces and moments at the 
boundary: 

(I-U b )0 b -{I-U a ) d a =T h a 

where a and b label the ends of the piece of beam, 
T h a are the resulting dead forces and moments 
applied in the interval [a,b\, and matrix U rep- 
resents the arm of the internal forces in the equi- 
librium equation of the moments: 

“ [ ~(P(0 ~t$x 0 

being x the pole the moments are referred to. 
The internal forces are defined by means of 
a constitutive law in terms of the generalized 
strains, i.e. d = d (ip) where d = {rpp} are the 
internal forces and moments and ip = {e,K.} are 
the strains and the curvatures. An arbitrary, 
complete beam section characterization can 
be used, that fully couples the deformations 
and the forces, thus allowing the modeling 
of anisotropic beams [6], [4], A three node 
parabolic C° beam element, that gives the 
exact solution for end- applied loads, has been 
implemented [15], [5]. The strains and curva- 
tures at the boundaries of the finite piece of 
beam are expressed as functions of the nodal 
configuration by means of shape functions, thus 
resulting in a finite element-like discrete beam. 
The finite volume beam can be regarded as 
a constraint that relates the reaction forces 
to the deformation of the link, and thus to 
the configuration of the system. Provided the 
relation between reactions and configuration 
is invertible (i.e. the Hessian of the strain 
energy is positive definite), the constraint 
equation can be implicitly solved, thus allowing 
the direct writing of the contribution of the 
beam to the equilibrium equations in terms of 
position and rotation unknowns. Finite volume 
beams are easy to implement in a multi-body 
formulation since only collocated evaluation of 
the contributions to the equilibrium equations 
is required. Moreover, they straightforwardly 
resemble the natural partition in distinct bodies 
that is peculiar to the multi-body formulation. 
The finite volume description of the deforma- 
tion of slender bodies is consistent with the 
mathematical, intrinsically discrete, multi-body 
model, and thus allows an easy but thorough 



modeling. 

Aerodynamic Forces. The aerodynamic 
forces are based on the strip theory, using el- 
ements that refer to rigid or beam shaped aero- 
dynamic surfaces. The aerodynamic coefficients 
are based on the interpolation of experimental 
data spanning 360 degrees of angle of attack. 
Corrections are made to determine the drag due 
to spanwise flow, as well as the effects of dynamic 
stall [7], Rotor elements are defined, to account 
for the effect of rotor induced velocity with an in- 
creasing degree of refinement, from uniform up 
to dynamic inflow modeling [17]. The imple- 
mented aerodynamic model is relatively poor, 
but it is satisfactory at least with regard to sta- 
bility analysis. The code is being interfaced to 
a wake analysis program for more accurate air- 
loads prediction. 

Tiltrotor Models 

The WRATS wind tunnel model is a semispan, 
1/5 scale aeroelastic model of the V-22 tiltro- 
tor [19]. It includes the right-side rotor, the 
pylon, and the half-wing. The rotor is pow- 
ered by a water-cooled electric engine; the con- 
trols are applied to the swashplate by means of 
three hydraulic actuators. The basic geometric 
and structural properties have been taken from 
Ref. [16], from the drawings and from direct 
measures on the wind tunnel model. Each sub- 
part of the tiltrotor has been modeled and anal- 
ysed in its basic kinematic and dynamic features, 
then the parts have been assembled together and 
the complete system has been analysed. By us- 
ing the same code and the same modeling for 
the single parts and the assembly, and by using 
a rather general approach in the kinematic and 
mechanical description of the single parts, any 
undue approximation has been avoided. The 
complete model is sketched in Figure 1. 

Submodels 

Blade. The single blade model has been 
used to analyse the dynamic properties of the 
blade, such as frequencies and aerodynamic 
properties. The flexbeam is described by a 
three-node beam element; a flexible control link 
is attached to the root of the blade by a rigid 
cuff and a pitch horn. The blade is joined to 
the outer end of the flexbeam by a spherical 
hinge, and by a spanwise oriented in-line joint, 


2.2 in. outwards from the rotor axis. These 
bearings allow the free pitching of the blade; 
the flap and lead-lag motions result from the 
bending the flexbeam, simply supported by 
the bearings. A rigid blade has been used for 
preliminary evaluation of the kinematics of the 
rotor; most of the analyses has been performed 
with a flexible blade model. Two discretisations 
have been used, made of two and four beam 
elements respectively. The flexibility of the 
flexbeam and the position of the joints heavily 
affect the stability of the blade. The flapping 
is dominated by the gimbal, whose motion 
determines a negative S 3 angle of about 15° (the 
pitch link is behind the blade). This is required 
to ensure the stability of the blade, being the 
rotor stiff-in-plane. But the virtual flap hinge 
due to the bending of the flexbeam is slightly 
outwards from the pitch link attachment, thus 
introducing a slight destabilizing effect on 
the flapping of the blade, which may become 
significant at high forward speeds, when the 
rotor acts as a propeller and the collective pitch 
at 75% of the blade is greater than 50°, because 
the flap and lag frequencies are very close to 
each other. This problem required a detailed 
study of the stiffness properties of the flexbeam. 

Gimbal. The hub is attached to the mast by a 
gimbal joint. This joint allows the hub to flap 
freely about an arbitrary axis normal to the 
mast, while the torque about the shaft axis is 
transmitted. The gimbal mechanism is made of 
two universal joints, respectively linked to the 
mast and to the hub at one arm of each cross. 
The other arms of the two crosses are connected 



Figure 1: Analytical Model 


to each other by a linkage, that transmits the 
torque between the mast and the hub and keeps 
the two universal joints at a constant distance. 
An in-line joint connects the hub to a spherical 
joint on the mast that represents the center of 
rotation of the hub. The gimbal is a complex 
joint that only in a very first approximation can 
be described by a simple universal joint. The 
multi-body approach easily allowed to model 
the complete kinematics of the joint from the 
beginning, thus ensuring the correctness of 
the formulation and with little computational 
overhead. The gimbal allows the rigid flapping 
of the whole rotor and, since the direction of 
the angular velocity tilts together with the hub, 
no Coriolis forces due to this motion result in 
the blades when the flapping is steady. At the 
same time, no stiffness due to centrifugal effects 
results from the 1 per rev. flapping motion, but 
that provided by a set of springs. 

Swash Plate. The blade pitch is controlled 
by a swashplate. It has been modeled by two 
rigid bodies that represent the plates, linked 
together by a plane revolution joint. The fixed 
plate is forced to slide along an in-line joint 
with respect to the pylon; its position and 
attitude are controlled by three actuators. The 
controls are imposed by changing the length of 
the actuators. The two plates are respectively 
linked to the helicopter and to the shaft by 
scissors, that constrain the rotation about the 
shaft axis. The control links of the blades 
are attached to the rotating plate by simple 
pins. This allows the exact modeling of the 
kinematics of the pitch controls. Extensive 
analyses of the resulting pitch-flap and pitch-lag 
couplings have been performed, and good 
correlation with available data and with di- 
rect measurements on the model has been found. 

Wing-Pylon. The wing model consists in two 
beam elements along the wing spar, at about 
25% of the chord from the leading edge; the con- 
version hinge joint is at the wing tip, at about 
75% of the chord, and it is elastically connected 
to the wing elements. The pylon is modeled as 
a rigid body. It is connected to the wing by 
the conversion joint, a flexible spindle that al- 
lows the free rotation about the conversion axis, 
while elastically constraining the bending in the 
other directions. The bending of the pylon about 
the axis normal to the wing is significant since it 
determines a coupled wing in-plane/pylon mode 


called the yaw mode, whose frequency is com- 
paratively low [18]. The conversion is controlled 
by a conversion actuator located ahead of the 
hinge; this actuator at present is not mounted on 
the wind tunnel model, and a downstop spring 
is used to emulate its stiffness properties. Dif- 
ferent sprins are used to change the conversion 
angle and to match the stiffness properties in 
the desired configuration. The analytical model 
has been tuned to reproduce these different stiff- 
ness configurations; the conversion actuator has 
been modeled, too, and conversion maneuvres 
have been simulated. 

Results 

The single submodels have been validated 
separately. Comparisons with experimental 
data and with results from previous analyses 
are presented. Partial assembly models of the 
rotor have been used for initial validations. 

Blade Models. The very first analyses 
consisted in checking the kinematic couplings 
of the blade. Sample results are reported in 
Figure 2 where a comparison is made between 
the computed and the measured pitch-flap 
coupling. Subsequent analyses involved the 
identification of the frequencies and of the 
modes of the clamped blade. A very basic 
comparison was made with a corresponding 
NASTRAN model of the blade, sharing the 
same properties and the same discretisation 
used for the multi-body analysis, to assess the 
correctness of the formulation of the beam 
element. The results are summarized in Table 1 

Rotor Models. The structural properties of 
the rotor have been analysed by means of many 



Collective Pitch at 75% of the Blade, deg 


Figure 2: Pitch-flap coupling 



Table 1: Cantilevered blade frequencies, Hz 


Table 3: Rot. freq., 888 rpm, # 7g % = —3°, Hz 


MODE 

GVT 

UMARC 

NASTRAN 

MBDyn 

4 el. 2 el. 

1 B 

12.29 

12.3 

11.5 

11.3 

11.7 

1 C 

34.11 

34.1 

33.4 

33.1 

32.7 

2 B 

52.44 

53.0 

56.7 

55.8 

55.0 

1 T 

113.35 

111.4 

127.0 

119.0 

122.0 


different partial models. Basically, data from 
Ground Vibration Tests (GVT) of the non- 
rotating assembly, with clamped pylon and free 
or locked gimbal, as well as computed results 
from various sources were available. The sin- 
gle blade with the flexbeam and the locked con- 
trol chain has been used to simulate the clamped 
gimbal condition. A complete model of the rotor 
has been used to simulate the free gimbal con- 
dition. It is important to notice that the three 
blade configuration and the gimbal joint break 
the symmetry of the system, thus requiring a 
modal analysis of the full rotor, being insuffi- 
cient a simple superimposition of the clamped 
rotor modes to the flapping motion allowed by 
the gimbal. The results from the multi-body 
analysis of the complete three-blade rotor were 
later confirmed by modifying UMARC to allow 
the modeling of the full rotor in the modal analy- 
sis phase, as detailed later. Results are reported 
in Table 2 

In vacuo rotating analyses have been compared 
only to numerical data since no experimental 
results were available. Myklestad and UMARC 
analysis results have been used. The full flight 
envelope has been tested, by performing modal 
analyses at the hover and forward flight rotation 
speeds and at different values of collective pitch. 
The change in pitch is significant when going 
from the hover condition to the forward flight 
at high speed, being of the order of 60°. Good 
agreement between all the available results has 


Table 2: Full rotor, non-rotating, Hz 


Mode 

GVT 

MBDyn 


10 deg. 

10 deg. 

50 deg. 

Gimbal 

2.0 

1.8 

1.5 

Cone 

6.8 

7.0 

7.8 

2 Flap 

25.0 

26.5 

36.5 

2 Flap asym. 

64.2 

57.1 

55.0 

3 Flap 

76.2 

78.0 

82.5 

1 Lag 

19.7 

19.0 

12.7 

2 Lag 

91.3 

98.0 

92.0 

1 Torsion 

112.1 

109.0 

107.5 


Mode 

Myklestad 

UMARC 

MBDyn 

Gimbal 

- 

14.8 

14.8 

Cone 

17.2 

17.3 

17.5 

1 Lag 

22.4 

20.8 

24.0 

Coll Lag 

42. 

44.0 

36.0 

2 Flap 

37.33 

49.6 

41.0 

2 Flap a. 

- 

70.2 

65.0 

3 Flap 

75.33 

90.3 

73.0 

Flap/Tors. 

89.33 

92.7 

90.0 

Lag/ Tors. 

- 

113.4 

104.0 

Torsion 

- 

116.0 

110.0 


Table 4: Rot. freq., 742 rpm, 0 75 % = 55°, Hz 


Mode 

Myklestad 

UMARC 

MBDyn 

Gimbal 

- 

12.4 

12.6 

Cone 

14.7 

14.9 

15.1 

1 Lag 

15.3 

15.8 

16.5 

2 Flap a. 

- 

42.3 

44.2 

Coll Lag 

32.7 

45.9 

46.9 

2 Flap 

45.3 

45.6 

49.1 

3 Flap a. 

- 

46.9 

60.3 

3 Flap 

66.0 

60.1 

65.2 

Flap/Tors. 

89.3 

90.6 

97.8 

Lag/ Tors. 

90.0 

90.8 

89.7 

3 Lag 

- 

92.0 

92.9 

Torsion 

- 

116.0 

108.5 


been found, as may be seen in Tables 3, 4. Data 
from Myklestad analysis refers to a clamped 
gimbal rotor, the missing values referring to 
non-symmetric modes that involve some gimbal 
flapping, as shown by multi-body and modified 
UMARC analyses. 

Wing Model. A relatively coarse mesh has 
been used for the wing, so only the very first 
modes have been sought, by comparing to GVT 
and to NASTRAN [16] detailed analyses of the 
model, see Table 5. Particular care has been 
taken in modeling the downstop springs that are 
used to simulate the stiffness of the conversion 
actuator in the different configurations. 

Wing-Rotor Models. The models of the wing 
and of the full rotor have been coupled to obtain 
the full detailed model of the tiltrotor. The com- 
plete model, with structural nodes, reaction un- 
knowns and control-related unknowns has 580- 
GOO degrees of freedom, depending on the type of 
analysis. A final check of the frequencies of the 
assembly has been made, followed by the aeroe- 
lastic analyses. The first wing modes are not 



directly affected by the modeling of the flexibil- 
ity of the rotor. The torsion mode of the wing 
is very close, and at some airstream speeds co- 
incident, to the rotor speed; this gives raise to 
resonance that can be seen in the frequency anal- 
yses of the internal forces of the wing. Four wing 
modes are mainly considered: the beam, chord 
and torsion modes of the wing, and the previ- 
ously mentioned pylon yaw mode. Referring to 
the fixed frame, the retreating rotor modes in- 
teract with the wing modes. This can be clearly 
appreciated from a frequency analysis of the 
wing response when the rotor modes are excited. 
Most of these modes cannot be easily identified 
when the aerodynamics are modeled, since they 
are highly damped. For this reason, a compre- 
hensive analysis of the structural properties of 
the model has been performed by simulating in 
vacuo operations, while the aeroelastic proper- 
ties have been estimated in different ways. The 
damping of the wing modes in forward flight 
has been estimated by system identification of 
the response to a given input, as is usually done 
during actual wind tunnel tests, while the aeroe- 
lastic pitch-flap coupling has been estimated by 
measuring the phase shift between an harmonic 
control input and the flapping response. Fig- 
ure 3 refers to a collective pitch maneuvre. It 
shows the geometric pitch of one blade as the 
collective control is raised from 0 to 10 degrees in 
one second. The simulation is performed in he- 
licopter mode; the nominal hover rotation speed 
of 888 rpm is reached in one second to obtain 
a trimmed condition (not shown); the residual 
oscillations at the beginning are due to the spin- 
up transient. There is no airstream speed. The 
difference between the given control and the ac- 
tual pitch of the blade is due to the deforma- 
tion of the flexbeam and of the flexible link. 
Figure 4 refers to the conversion maneuvre per- 
formed by a deformable blade model. It shows 
the internal forces at the wing root. The conver- 


Table 5: Wing frequencies, Hz 


Mode 

d/s 

GVT 

NASTRAN 

MBDyn 

Beam 

on 

6.00 

6.16 

5.9 


off 

5.51 

5.45 

5.4 

Chord 

on 

8.45 

9.33 

9.1 


off 

8.45 

8.74 

00 

bo 

Torsion 

on 

12.5 

12.6 

12.5 


off 

10.6 

10.6 

11.0 

P. Yaw 

on 

16.5 

18.9 

17.2 


off 

16.7 

16.7 

16.6 


sion is performed at a 10 deg/s constant angular 
speed. Oscillations of the internal forces due to 
the untrimmed initial conditions are appreciably 
damped as the maneuvre proceeds to the end, at 
t = 9 s. The following abrupt raise of the inter- 
nal forces is due to the transient caused by the 
sharp end of the maneuvre. The flexible blade 
model has been used to simulate the response 
to a cosinusoidal vertical gust in airplane mode, 
of 10 ft/s amplitude. Both the stability and the 
sensitivity of the tiltrotor have been addressed. 
Figures 5, 6 show the wing out-of-plane internal 
moment due to the gust at different airstream 
speeds, for both off- and on-downstop configu- 
rations. In figure 5 the off-downstop configura- 
tion is clearly less damped than the other one, 
in fact the stability boundary in air is about 160 
Kts. When the rotating speed is increased, the 
stability boundary moves towards lower speeds, 
as shown by previous analyses and experiments 
[18]- 

Computational Notes 

The complete model has nearly 600 degrees of 
freedom. No modal condensation has been per- 
formed; the physical flexible elements have been 
used throughout the analyses. A typical model 
of the tiltrotor is made of 45 nodes, 39 rigid 
bodies, 35 joints of different kind, 18 beam el- 
ements, 14 aerodynamic elements, 6 control- 
related nodes and 4 control-related elements. 
The simulations have been performed on off- 
the-shelf PCs. The time step initially required 
for the rigid blade model was At = .5 x 10 -3 
s, while the deformable blade model required 
At = .25 x 10 -3 s to start correctly. When 
a variable time step was used, the rigid blade 
model simulations quickly reached a value of 



blade pitch 

control input 



0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

t, sec 

Figure 3: Collective pitch maneuvre 




3.0 3.5 x 10 -3 s, while, the deformable blade 

model ones, reached about 1.0 1.2 x 10 -3 s. 

The conversion simulation required about 4.5 
hours on a Pentium PRO 200 for a total of 40000 
fixed size time steps (10 s at At = .25 x 10 -3 
s). When performed with variable step size, it 
required about one hour. After the model was 
refined, and a soft start was used, the flexible 
blade model is able to start with At = 10 -3 
s, requiring about 1.7 hours, or 1.1 hours on a 
Pentium II 333 and 0.8 hours on a Pentium II 
450. Tests are being performed on Digital work- 
stations with Alpha processor. The speed has 
been increased of a factor 4.5 for typical simula- 
tions. These numbers make this kind of analy- 
sis interesting even for a large, time consuming 
parametric study for rotorcraft design, not only 
for analysis. 

Concluding Remarks 

The paper illustrates the feasibility of a multi- 
body, global modelling approach to the analy- 
sis of a rotorcraft. A tiltrotor has been chosen 
since its peculiarities are likely to highlight the 
possible limitations of more conventional anal- 
ysis and design formulations when facing non- 
conventional configurations that may require the 
analysis of non-conventional flight conditions or 
maneuvres. The formulation proved to be effi- 
cient without excessive simplifications. The effi- 
ciency has been preserved while maintaining the 
physical meaning of both the equations and the 
unknowns. The multi-body formulation relies on 
the writing of the equilibrium equations of each 
body of the system, regarded as independent. 
Kinematic constraint equations are added, re- 



Figure 4: Conversion — Bending moment 



Figure 5: Gust — off-downstop 


suiting in reaction unknowns that contribute to 
the equilibrium of the bodies. Finite, large ro- 
tations are accounted for in a computationally 
efficient way, based on an updated Lagrangian 
approach, that assumes the predicted configu- 
ration of a node as reference. The writing of 
the Jacobian matrix is simplified, resulting in 
computational time saving, while the accuracy 
is preserved by consistently computing the resid- 
ual. A model of the tiltrotor used in WRATS in- 
vestigation has been analysed, consisting in ro- 
tor models of increasing refinement, with rigid 
and flexible blades, the gimballed constant ve- 
locity joint, the swashplate and the control links, 
the pylon, the conversion hinge and the flexible 
wing. Aerodynamic loads have been considered, 
to simulate different test conditions, from aeroe- 
lastic stability investigations to the simulation of 
complex manoeuvers, including conversion and 
blade pitch control. The results here presented 
should be considered as an assessment of the 
validity of the formulation, rather than a com- 
plete analysis of the WRATS tiltrotor model. 
The multi-body code that resulted form the pro- 



Figure 6: Gust — on-downstop 






posed formulation is currently being improved 
by adding true active control capabilities, inte- 
grated electric and hydraulic networks analysis, 
modal description of flexible bodies, and, from a 
computational point of view, iterative, matrix- 
free solvers to speed-up the solution phase in 
view of a low-cost, coarse scale parallelization of 
the code. 
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